Magnetometry with nitrogen-vacancy ensembles in diamond based on infrared 
absorption in a doubly resonant optical cavity 
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We propose to use an optical cavity to enhance the sensitivity of magnetometers relying on the 
detection of the spin state of high-density nitrogen-vacancy ensembles in diamond using infrared 
optical absorption. The role of the cavity is to obtain a contrast in the absorption-detected magnetic 
resonance approaching unity at room temperature. We project an increase in the photon shot- 
noise limited sensitivity of two orders of magnitude in comparison with a single-pass approach. 
Optical losses can limit the enhancement to one order of magnitude which could still enable room 
temperature operation. Finally, the optical cavity also allows to use smaller pumping power when 
it is designed to be resonant at both the pump and the signal wavelength. 

PACS numbers: 76. 30. Mi, 78. 30. Am, 07.55. Ge, 42.60.Da 



I. INTRODUCTION 

The negatively charged nitrogen- vacancy (NV~) cen- 
ter in diamond can be used as a solid-state mag- 
netic sensor due to its electron spin resonance (ESR). 
The center can be optically polarized and its polariza- 
tion detected through the spin-state dependence of the 
luminescence^. Sensors based on a single NV~ cen- 
ter have the potential to achieve atomic-scale spatial 
resolution 3 5 . On the other hand, magnetic field sensitiv- 
ity can be enhanced by engineering the diamond material 
in order to increase the spin dephasing time which limits 
the ESR linewidth 6 .. The magnetic response of an ensem- 
ble of NV~ centers 7 - - — leads to a luminescence magnified 
by the number N of the sensing spins. Such collective 
response also improves the signal to noise ratio and the 
sensitivity by a factor \/N since the quantum projection 
noise associated with the spin-state determination scales 



Currently, the sensitivity of practical magnetometers 
based on the detection of red luminescence of the NV~ 
ensemble is limited by background fluorescence and 
poor collection efficiency. Recent advances in diamond 
engineering have enabled improvements in collection 
efficiency which should improve fluorescence based 



sensors 



12-16 



but here we consider a different approach. 
In addition to the well-known transitions leading to red 
fluorescence, it has been shown recently the existence 
of an infrared (IR) transition related to the singlet 
state s 17 ! 18 . This transition can be exploited in an 
IR-absorption scheme with an increased sensitivity as 
compared to the usual scheme^. In this paper we show 
that using IR absorption detection in combination with 
a high-finesse optical cavity, it is possible to tune the 
absorption contrast to order unity thereby dramatically 



improving the magnetic field sensitivity. We first recall 
the parameters which set the magnetometer sensitivity. 
We then theoretically investigate the extension of 
this detection scheme to the case where the diamond 
crystal hosting the NV~ ensemble is inserted inside a 
high-finesse optical cavity, as it is usually done in cavity 
ring-down spectroscopy 20 . Finally we determine the 
improvement of the magnetometer response associated 
with the cavity quality (Q) factor. 



II. SINGLE-PASS PHOTON SHOT-NOISE 
LIMITED MAGNETIC FIELD SENSITIVITY 

The principle of the method is similar to the one used 
in optical magnetometers based on the precession of spin- 
polarized atomic gases 21 . The applied magnetic field 
value is obtained by optically measuring the Zeeman 
shifts of the NV~ defect spin sublevels via the absorption 
monitoring of the IR probe signal. The photodynam- 
ics of NV~ centers are modeled using the level structure 
depicted in Fig [TJa). The spin sublevels m s = and 
m s = ±1 of the 3 A2 ground triplet state are labeled |1) 
and 1 2) and separated by D = 2.87 GHz in zero mag- 
netic field. 1 3) and |4) are the respective spin sublevels 
of the 3 E excited level. Levels |5) and |6) are single- 
state levels related to the infrared absorption transition. 
The relaxation rate from state i to j is denoted kij. As 
&35 <C (see Table [J in Appendix [Ej), the system is 
optically polarized in m s = while pumping the NV~ 
centers via the phonon sideband. Without microwaves 
applied, there is reduced population in the metastable 
singlet state, |6), corresponding to a minimal IR absorp- 
tion signal. Under application of resonant microwaves 
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with frequency D ± jB/(27r), where B is the magnetic 
field projection along one of the four NV - orientations 
and 7 = 1.761 x 10 11 s _1 T _1 is the gyromagnetic ratio, 
population is transfered from m s = to m s = ±1 sub- 
level resulting in greater population in the metastable 
singlet and lower IR signal transmission. The experi- 
mental configuration for single-pass absorption measure- 
ments is shown in Fig[TJb). The output transmission is 
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FIG. 1. a) Level structure of NV - center in diamond. The 
photophysical parameters related to this six-level system are 
given m Tab. [T]of Appendix [H The solid (dot) lines corre- 
spond to radiative (non-radiative) transitions. D « 2.87 GHz 
is the zero-field splitting of the ground state, b) Diagram of 
the experimental configuration used to measure the single- 
pass contrast of the IR absorption under resonant microwave 
application^. 7o,p and io,s are the pump (wavelength Ap) 
and the probe input intensities. 

measured either with or without applying the resonant 
microwaves. The contrast C is defined as the relative dif- 
ference in the IR signal detected after propagation in the 
diamond crystal of length L 



C = 



Iout,s(0) - [out, s(^r) 
/out, 5(0) 



(1) 



where / O ut,s(0) [lout, s(^r)] denotes the IR signal inten- 
sity without [with] the application of the microwave field 
whose Rabi angular frequency is denoted Qr. We can es- 
timate the photon shot-noise limited sensitivity at room 
temperature for an optical power compatible with the 
IR saturation intensity. For an ESR full-width-at-half- 
maximum (FWHM) r mw , the magnetic field sensitivity 
(or the minimum detectable magnetic field) of a mag- 
netometer based on IR absorption measurement is given 
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where Ps is the measured IR probe beam signal out- 
put power (wavelength As), and t m is the measurement 
time. Assuming no power broadening from either pump 
or microwaves, the ESR FWHM is related to the electron 
spin dephasing time by r mw = 2/T 2 * (in rad/s). For a 
detected IR signal power Ps = 300 mW using Eq. (|2j) 
with parameter values given in Tab. [TT] of Appendix [E] 
we obtain a shot-noise limited magnetic field sensitiv- 
ity of 20 pT/^/Hz in a single-pass configuration at room 
temperature. Note that considering this IR signal power 
and a beam waist diameter of 2wo = 50 fim there is 
no saturation of the IR absorption (see Appendix |B|) . 
For this single-pass configuration, the contrast cannot be 
improved by increasing the thickness of the sample since 
for L larger than the pump penetration depth (« 120 jam 
from the absorption cross section and NV - center density 
of Tab. HI and HI]) its absorption becomes too strong. The 
photon shot-noise limited sensitivity can be compared to 
the spin-noise limited sensitivity 



SB q = 



(3) 



where we take into account through the factor of 2 that 
only one fourth of the NV - centers are oriented along the 



magnetic field 



is the NV -center density and V is 



the illuminated diamond volume. In the single pass con- 
figuration of Ref. [19], the spin- noise limited sensitivity 
is about 0.02 pT/VHz. 



III. SENSITIVITY ENHANCEMENT 

According to Eq. ((2]), the magnetic-field sensitivity is 
limited by the low contrast C. In particular, at room tem- 
perature the contrast is an order of magnitude smaller 
than at 75 K due to homogeneous broadening 19 . It can 
also be seen as limited by the optical depth estimated 
to only 2.2 x 10 -2 for the experimental demonstration 
reported in Ref. [19]. However, the optical depth can 
be increased by using a cavity resonant at the IR signal 
wavelength resulting in an increase of the optical path by 
a factor proportional to the finesse of the cavity. More- 
over, using a diamond crystal thickness smaller than the 
pump absorption length allows to overcome the issue of 
the pump depletion and to obtain a good microwave field 
homogeneity along the crystal. We consider the Fabry- 
Perot cavity configuration depicted in Fig [2^i), consist- 
ing of a two-side coated bulk-diamond plate containing a 
high NV - -center density (larger than 4 x 10 23 m -3 ). We 
consider an all-pass Fabry-Perot cavity for the IR signal. 
This means that the amplitude reflectivity of the back 
mirror is Pback,^ = 1 and of the input mirror reflectivity 
is Pin, s < 1- Regarding the pump, we consider either 
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FIG. 2. a) All pass cavity (we consider a perfectly reflecting 
backside mirror |pback,s| = 1) used for magnetic field sensi- 
tivity enhancement. pi n ,i is the amplitude reflectivity of the 
input coated mirror. The cavity can be doubly resonant for 
the pump and the signal. C: optical circulator, F: optical 
filter rejecting the pump beam, D: optical detector, b) Re- 
flected spectrum from the cavity for switched-on or switched- 
off microwaves (mw) resonant at the level |1)-|2) transition. 
An is the IR cavity resonance wavelength, c) Maximal intra- 
cavity pump beam optical power magnification factor for a 
given value of intracavity absorption and two values of cavity 
lengths L — Lo and L — 10Ln. The finesse of the cavity at the 
pump wavelength is denoted Fp. Lo is the cavity length which 
gives a critical coupling (and thus the the optimal magnifica- 
tion factor) for Fp « 100. Note that even with pi n ,p — 0, the 
all-pass configuration gives a maximal magnification around 
4 due to reflection on the backside mirror. 



single-pass propagation (pi n ,p — Pback,p = 0) or all-pass 
cavities (pback,p = 1)- We define the reflection of the 
cavity at optical resonance by R{ = I out with (with 
ie{P,S}). 



pendix [Cl at the first order 
Rs(0) = 
Rs(n R ) = 



a 



a 



a n R 



(4) 



The finesse of the cavity given in Eq. (|C5|) can also be 
written at the first order in e and a* 



(5) 



where i — for off-resonance microwaves and i — VLr for 
on-resonance microwaves. 



1. Optimal cavity coupling 

Assuming a perfect spin polarization and no additional 
optical losses, we have: ao = 0. In this case, Rs(0) = 
1 and thus the off-resonance reflected detected signal is 
equal to the input signal power Po,5- The contrast reads 
C = 1 — Rs(Qr) and the magnetic field sensitivity is given 

by 
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For s = d£i R1 the incoming and outgoing fields destruc- 
tively interfere at the resonant wavelength and Rs(^r) — 
0. The laser probe beam is then critically coupled 25 
to the cavity-NV _ ensemble system and the contrast is 
equal to 1 . For this particular value the optimal sensibil- 
ity of the magnetometer is reached. 



2. Effects of the microwave off-resonance absorption 



A. Basic principle of the cavity effect 

The complete analysis of the cavity has to be per- 
formed numerically. In order to allow a simple inter- 
pretation of the results, we first derive analytical expres- 
sions for the sensitivity assuming no saturation of the 
IR-signal absorption. The absorption of the IR signal 
due to levels |5) and |6) and the spin polarization due 
to the pump beam system is simply taken into account 
by As the single-pass round-trip amplitude transmission. 
We also assume a good finesse cavity at the IR signal 
wavelength and thus the input mirror reflectivity can be 
written Pin,^ = 1 — s with e <C 1. With the application of 
the resonant microwave field we have: As(£Ir) = l — an R 
( a Q R <^ 1) whereas for an off-resonance microwave field 
we have: As(0) = 1 — ao (ao <C 1). We define the op- 
tically resonant reflectivity for respectively off- and on- 
resonance microwave fields using the results given in Ap- 



Now we consider the more realistic case of a non ideal 
spin-polarization and material with parasitic IR losses 
which gives clq r > ao > 0. There are three possible cases 

i) e > ^/a~o~a~^ R s (0) > R S (^r) 
< ii) £ = v/«^ Rs(P) = R s (n R ) (7) 
k iii) e < ^a an R Rs(0) < R s (£Ir). 

Consequently, depending on the relative value of Rs(0) 
and Rs(Qr), the expression of the contrast is different. 
This can be taken into account by writing 

= \R S (0) - Rs(n R )\ ( ) 

m &x [Rs(n R ),Rs(0)Y { ' 

This relation can be used to write the expression of the 
minimum detectable magnetic field taking into account 
the detrimental effect of the residual IR absorption due 
to non-ideal branching ratio to the metastable state by 
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multiplying P ^ by max [#5 (Or), i2s(0)] to obtain the 
detected IR power Ps of Eq. ([2]). The fundamental ad- 
vantage of the present method is that this quantity falls 
under the square root whereas for methods based on the 
visible-fluorescence monitoring the non-ideal branching 
ratio reduces the contrast C by a similar amount, but 
this quantity falls outside the square root. One can esti- 
mate that in the same conditions, the minimal detectable 
magnetic field SB is reduced by a factor of « 5 in com- 
parison with SB? obtained via fluorescence method with 
a collection efficiency 77 « 0.47 (see details and discussion 
in Appendix IF]) . The sensitivity thus reads 



SB = 



I hex max [Rs(n R ),Rs(0)} 



j\Rs(0)-R s (n R )\ 



Po,st m \s 



(9) 

In the present case, there are two critical-coupling condi- 
tions, thus the sensitivity 5B can reach two optimal val- 
ues obtained for e = a^ R (solid line in FigEfc) or e = clq. 
Note that due to the factor y^max ^(0)] in 

the numerator of Eq. (|9]), the minimum values of 5B is 
actually reached for values of e slightly different from the 
exact critical-coupling finesse. This will be accurately de- 
scribed in the numerical calculations. We first consider 
the case i) of Eqs. (0). Assuming e ^> an R we have 



SB 



7rT n 
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4ryFs(an R - ao) V Po,st m ^s ' 



(10) 



This means that for low cavity finesses the effect of the 
cavity is to reduce the minimum detectable magnetic field 
value by a factor equal to the finesse Fs ~ tt/s. For 
£ = v" a ° a ^ (case ii), the contrast is equal to zero and 
SB reaches a singular value as shown in Eq. (|2]) . Finally, 
for e < ^aoan R (case iii), assuming e <C ao the sensitivity 
reads 
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(11) 



This shows that the sensitivity can be greatly impaired 
(i.e. SB increases) if the empty cavity finesse (ir/e) is 
larger than that of a critically coupled cavity given by 
7r/(2a ). Moreover, Eqs. ([TOj) and ((TTJ) show that if the 
off- and on-resonance loss values ao and clq r are too close, 
the sensitivity is also impaired. 

As a conclusion, the level |6) is always partly popu- 
lated due to the non ideal branching ratio to the dark 
singlet state (ks5 7^ 0). This results in absorption of 
the IR probe beam, even in the microwave-off state (i.e. 
no resonant microwaves applied) and the implementa- 
tion of a cavity will also increase this effect and reduce 
the detected IR photon number I ou t,s- Thus, the cav- 
ity induces simultaneously an increase in the contrast C 
and a reduction of the detected photon number in the IR 
beam. Consequently, for a given single-pass absorption, 
the cavity finesse cannot be arbitrarily increased and the 
magnetic field sensitivity SB reaches a minimum value 



intrinsically limited by NV _ photophysical parameters 
and by diamond intrinsic IR optical losses. Those effects 
are quantitatively described in the next section where 
numerical results are reported. 



B. Numerical calculations 

The output fields E out ^ both for the pump and IR 
signal are deduced from the input and intracavity for- 
ward and backward propagating fields f%(z) and bi(z) 
described Fig. [2ji) using the slowly varying envelope ap- 
proximation. Note that the intracavity absorption (ob- 
tained by solving the six-level rate equations) depends 
nonlinearly on the intracavity intensity U(z) and thus a 
numerical optimization routine on fi(L) must be used 
to deduce the reflected powers both at pump and signal 
wavelengths for the target values 26 of io,i (see Appendix 
Ofor details on the calculation method). 

We consider two NV~ center concentrations^ 2 ^ i) con- 
figuration 1: n = 4.4 x 10 23 m -3 and T 2 * = 390 ns ii) con- 
figuration 2: n = 28 x 10 23 m -3 and T 2 * = 150 ns. For 
high NV~ -center density, single-pass absorption is high 
and the system is less sensitive to parasitic optical losses, 
but the electron spin dephasing time is shorter than for 
less low density samples. For each of these configurations 
we analyze: i) the effect of the diamond crystal sample 
thickness, ii) the effect of the input power, and iii) that 
of the Q-factor of the cavity. The Q-factors are defined 
by Qi = 2n d LFi/\i, (i e {P, S}) n d = 2.4 being the dia- 
mond refractive index and where we recall (see Eq. (|C5|) 
in Appendix [C]) that the finesse Fi is defined by 



Fi 



1 - Pin,iAi 



(12) 



with Ai the single-pass round-trip transmission. Note 
that in the case of a resonant pump field, the cavity is 
designed in order to reach exactly the critical coupling 
Ap = Pin,p which gives the maximal intracavity pump 
field enhancement and the optimal pump energy transfer 
to the NV~ ensemble. 

Figure [3] shows the magnetic-field sensitivity as a func- 
tion of the cavity Q-factor Qs at the IR-signal wavelength 
for two cavity lengths and three values of as which repre- 
sents the IR-signal optical-loss due to the bulk diamond 
material alone. In the rate-equation approximation, the 
sensitivity reaches two maxima (minima of SB), the first 
corresponding to a cavity critically coupled when the mi- 
crowaves are switched-on and the second corresponding 
to a cavity critically coupled when the microwaves are 
switched-off. Between these two optimal coupling con- 
figurations, we observe a sharp decrease of the sensitiv- 
ity corresponding to a cancellation of the contrast. For 
this particular situation, the reflection for the microwave 
switched-on and switched-off cases are equal. The IR 
optical losses reduce the sensitivity of the cavity but for 
as = 0.5 cm -1 (as = 0.1 cm -1 ) the best sensitivity 



5 



Config. 1 - /7=4.4xl0 23 nrf 3 



100 ■ 
10 • 



1 1 — I — I I I I 1 1 1 1 — I — I I I I 1 1 

a) / =100|im / F =23 / / =80MW/m 2 



a =0.1cm 

s 

-a s =0.5cm 
= a =3cm 1 




10 • 



°o o.l- 



1 1 1 




■1 




— — 






l 


^^^^ J 






l 


b) /.=10nm / F p =160 / / p =8MW/m 2 









10' 



1(T 



io J 



Config. 2 - n=28xl0 23 m" 3 




' 1 2 

CQ : d) L=10|im / F =31 / / =40MW/m 
to n 1 " p p 

U.l — | 1 1 1 — i l l l l | 1 1 — l — l l l l l | 



10' 



10" 



10 



io 5 



Cavity quality factor Q s at ^ s 



FIG. 3. Shot noise limited magnetic field sensitivity vs Q- 
factor of the cavity at the signal wavelength and for dif- 
ferent values of IR-signal optical losses (as)- Calculations 
are done for tt R = 2tt x 1.5 MHz, P ,s = 300 mW with 
Jo,s = 150 MW/m 2 and no optical losses for the pump 
(ap = 0). c) Config. 2 and L — 100 /mi, we assume a single 
pass pumping. For each plot, the value of 5B obtained for low 
Q s is about half compared to that obtained for single-pass 
propagation as expected from the use of a high reflectivity 
backside mirror. 



can reach 0.6 pT/vHz (0.3 pT/vHz) corresponding to 
almost two orders of magnitude enhancement in compar- 
ison to single-pass approaches. For strong optical losses 
(as = 3 cm -1 ) the sensitivity is still enhanced by more 
than one order of magnitude and the performance of the 
cavity system is comparable with that of the same sam- 
ple in a single-pass configuration at low temperature^. 

We now discuss the results for IR optical losses set to 
as = 0.5 cm -1 . For n = 4.4 x 10 23 m -3 , it is possible to 
use a doubly resonant cavity to increase the intracavity 
optical pump intensity and thus to reduce the required 
external intensity as illustrated in Fig. [2k). By diminish- 
ing the length of the cavity, the single pass attenuation 
is reduced and thus it is possible to increase the pump 
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FIG. 4. Shot-noise limited magnetic-field sensitivity calcu- 
lated for Qr = 2?r X 1.5 MHz, as = 0.5 cm -1 , ap — and 
2w o = 50 /mi varying the input IR signal power. For Config. 2 
(n — 28 x 10 23 m -3 ) and L — 100 /mi we assume single-pass 
propagation for the pump. The cavity parameters have been 
optimized using Fig. [3] 



cavity finesse and thus to strongly reduce the required 
amount of pump power from 400 MW/m 2 (single-pass 
propagation) to 8 MW/m 2 . For n = 28 x 10 23 m~ 3 , the 
pump absorption is so high that for L = 100 jam a dou- 
bly resonant approach does not give any improvement in 
the required pump power (Iq,p = 400 MW/m 2 ). Never- 
theless, for short cavities (L = 10 fim) a modest-finesse 
cavity for the pump (Fp = 31) leads to a reduction of 
the external pump power (down to Jo, p = 40 MW/m 2 ). 
In Fig. [Hwe plot the magnetic field sensitivity as a func- 
tion of the IR signal input power Pq^s for a beam- waist 
diameter 2wo = 50 fim. For thick diamond slabs, the 
saturation is obtained at high power (> 10 W). For thin 
diamond slabs, the use of high-finesse cavities reduces 
the signal saturation power. In the highest- Q-factor case 
(Config. 1 and L = 10 /im), saturation starts around 
Po,s ~ 300 mW. For high signal input power thermal 
effects must be taken into account. Note that these ef- 
fects would improve the sensitivity via the thermo-optic 
effects. More generally any nonlinear dispersive effect 
would increase the sensitivity of the device. In this case, 
a change in the absorption for the signal would induce a 
shift of the cavity resonance. In the example of Fig. EJb), 
if we denote Ai — Ao the shift of the cavity, the contrast 
would be given by [R s (0,\i) - R s (tt R )}/R s (0, Xi) and 
would have approximately the same value than without 
nonlinear effects. However the detected reflected power 
would be Rs(0, Ai) x Pq^s and would be greatly increased 
in comparison with Rs(0) x P^s which could reduce the 
value of the minimum of the detectable magnetic field as 
shown for example by Eq. ([2]). 

We can check that all the results given here are consis- 
tent with the quantum-noise limited sensitivity: i) Con- 
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fig. 1 SB^ = 0.2 pT/vIS and SB^ = 0.06 pT/V^ ii) 
Config 2. 5B q = 0.13 pT/V^ and SB q = 0.04 pT/x/Fte 
for L = 10 /im and L = 100 /im respectively. The choice 
of parameters for each case considered above results from 
an optimization depending on the crystal thickness and 
NV~ center concentration. Note that in the most res- 
onant configuration (Config. 1 and L = 10 /im), the 
optimal overall Q-factor of the cavity for the probe is 
around 5.3 x 10 4 , giving a cavity bandwidth 7 cav = 
2tt x 5.4 GHz much larger than the probe-laser linewidth 
(jl ~ 27r x 10 MHz) used for single-pass experiments re- 
ported in Ref. [19j- For high NV~ concentrations (Con- 
fig. 2), the required Q-factor can be low (< 3 x 10 4 ) and 
thus the total optical path £ = XsQs / {^ n d) ~ 2 mm) 
is smaller or almost equal to the Rayleigh range obtained 
for a waist diameter 2wq = 50 /im (2Zr « 3.8 mm). Con- 
sequently, the simple planar Fabry-Perot geometry 28 de- 
picted in FigEk) can be used. Finally, considering highly 
concentrated thin samples the required Q-factor can be 
around 2 x 10 4 which is compatible with recent measure- 
ment reported on integrated diamond microcavities 29 . 

C. External-mirror cavities 

For the highest-finesse cavities, appropriate for a con- 
centration of n = 4.4 x 10 23 m -3 , the effective length £ is 
longer than the Rayleigh range for the chosen beam waist 
value (2wq = 50 /im). Consequently, external spheri- 
cal mirrors should be used. If we consider for example 
a confocal cavity, the distance between the mirrors is 
L cav = 2Zr — 3.8 mm. For a 100 /im (10 /im) thick dia- 
mond plate, the finesse of the cavity would be Fs = 110 
(Fs = 1150). Consequently, in the case of the highest fi- 
nesse cavity, the Q-factor would be 8.4 x 10 6 correspond- 
ing to a cavity bandwidth 7 cav = 2tt x 34 MHz still larger 
than the probe-laser linewidth. We have assumed here 
distributed optical losses such as as = 0.5 cm -1 ; if we 
consider that optical losses mainly come from diamond 
interface roughness, it implies that in the more unfavor- 
able case (for the 10 /im-thick diamond plate), the root 
mean square deviation of the surface to planarity of the 
diamond interfaces 3 -^ has to be less than 2 nm, which is 
attainable with state-of-the-art fabrication techniques 31 . 

IV. CONCLUSION 

The use of a cavity can enhance the sensitivity of opti- 
cal magnetometers based on IR absorption of NV _ cen- 



ters in diamond at room temperature. We found that for 
diamond samples with a high density of defects (NV~- 
center concentration larger than n > 4.4 x 10 23 m -3 ), 
our configuration allows an enhancement of two orders 
of magnitude in comparison with single-pass configura- 
tions. In the presence of high IR optical losses the en- 
hancement is reduced to one order of magnitude. The 
use of a cavity compensates for the reduction of the 
optical depth due to homogeneous broadening at room 
temperature^. Moreover, doubly resonant (for the pump 
and the probe) cavities can be used to reduce the amount 
of required pump intensity (down to 8 MW/m 2 ). Us- 
ing diamond samples with a very high density of de- 
fects (n ~ 28 x 10 23 m -3 ), this approach could be im- 
plemented using monolithic planar Fabry-Perot cavities 
or integrated diamond photonic structures such as mi- 
crodisk or microring resonators. For smaller defect con- 
centrations (n ~ 4.4 x 10 23 m -3 ), external spherical- 
mirror cavities should be used. 



ACKNOWLEDGMENTS 

This work was supported by the France-Berkeley Foun- 
dation, the AFOSR/DARPA QuASAR program, IMOD 
and the NATO Science for Peace program. K.J. was sup- 
ported by the Danish Council for Independent Research 
Natural Sciences. 



Appendix A: NV six-level modeling 

The local density rij(z) (with j G [1,6]) of the centers 
of each level are calculated by solving the rate equations 
assuming drtj/dt — 0. We consider spin-conserving op- 
tical transitions. The pump excites a vibronic sideband 
which decays quickly via phonon emission to levels |3) 
and |4). This allows us to neglect the down-transition 
rates due to the pump light. At z, the relation between 
the optical intensity and the center densities is given by 

M(z)-M{z)=M Q , (Al) 

where A/o = (0, 0, 0, 0, 0, n) T , M contains the values of 
the center densities: M = (ni, 77,2, ^3, ^4, ^5, tlq) t and 
the matrix A4(z) can be written: 



M(z) = 



[W P (z) + W mw ] 


w mw 


hi 








hi 


\ 


w mw 


-[W P (z) + W mw ] 










h2 


W P (z) 





-(*31 + ^35) 
















W P {z) 





-(*42 + ^45) 
















&35 


&45 


~[W s (z) 


f r] +w s (z) 




1 


1 


1 


1 


1 


1 


/ 
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we assume here a closed system: Y^j=i n j ~ n - The 
transition rates Wi (i = P for the pump and i = S 
for the IR signal) are related to the optical intensity 1$, 
the wavelength A$ and the absorption cross section ai 
by Wi = GiliXij (he). Assuming a low Rabi angular fre- 
quency Or, in the rate-equation approximation, the mi- 
crowave transition rate is calculated as W mw = ^T 2 */2 
where T 2 * is the electron spin dephasing time. The center 
density in each level is calculated by J\f = A^ _1 A/q. 



Appendix C: Analytic expression of the 
cavity-reflectivity in the linear regime 

Here we consider the cavity described in Fig. [2^) with 
Pback,^ = 1- We denote the probe input field ^0,5, the 
reflected field i^out,^ and the forward propagating field 
inside the cavity at the input mirror J~s(0). Introduc- 
ing the amplitude mirror IR transmission coefficient Ki n ,s 
verifying ft 2 n s + pf n s = 1 and the round-trip phase (p, 
we can write 



Appendix B: IR-absorption cross section estimation 

In order to model the system we have to evaluate the 
IR absorption (due to singlet states) cross section as 
which has not been measured so far. With the aim of de- 
signing a cavity based magnetometer, the value of as is 
important to evaluate the intracavity IR signal intensity 
saturation. This completes the already reported list of 
photophysical properties of the NV~ centers in diamonds 
that are summarized in Tab. [J given in Appendix [EJ Here 
we estimate as by using the single-pass IR-absorption 
measurements described in Ref. [19[. We assume that 
the measured magnetic field is oriented in such a way 
that the microwaves are only resonant with NV~ centers 
of a particular orientation, i.e., one quarter of all the 
NV~ centers 2 ^. In the single-pass configuration C can be 
calculated by integrating the two differential equations 
considering off-resonance pumping and a resonant exci- 
tation (including stimulated emission) for the signal 



dip 

dz 

dz 



{a P [ni(z) +712(2)] + ap}Ip(z) 
{a s [n 6 (z) -715(2)] +as}Is(z), 



(Bl) 



where the densities rii(z) with i G [1,6] are the station- 
ary solutions of the rate equations corresponding to Fig. 
□Ji) (see Appendix IX]) . with i G {P, S} are the op- 
tical losses due to light scattering or parasitic absorp- 
tion. Calculations are carried out using the parameters 
given in Ref. [HI recalled in Tab. |TT] (see Appendix 
IE]) . The two unknown values are the IR absorption cross 
section as and the optical losses ap at the pump wave- 
length. The method consists in numerically finding the 
values of as which gives the contrast value defined in 
Eq. (pQ) and reported in Ref. [l9|. We have then de- 
duced that for a monochromatic excitation (the linewidth 
of the IR laser is j L w 2tt x 10 MHz <C 7ir), the IR 
absorption cross section due to the metastable level is 
as = (2.0 ± 0.3) x 10 -22 m 2 . The uncertainties come 
from the value of ap which has been assumed to vary 
from to 10 cm -1 . The associated saturation intensity 
is J sa t,5 = hcT/(2X s a s ) « 500 GW/m 2 . 



^s(0) = j^sEo^s + Pin,sAsFs(0)e jip 
E ut,s = Pin,sEo,s +3K>m i sA s J r s(0)e J<p . 



(CI) 



By eliminating .Fs(O), we can deduce the amplitude 
transfer function of the cavity 



E ( 



out, S 



A s e j 



(C2) 



E ,S 1 - PmMs&V 

The intensity reflectivity of the cavity is thus given by 



E ( 



out, S 



Eq,s 



Pi 2 n,S 



A% - 2p^ s A s cos Lp 



1 + pfn,S A S ~ 2 Pin,sA S COS (p ' 



(C3) 



At resonance cp — (2ir) the reflectivity of the cavity can 
be written 



Rs = 



Pin,S 



1 - Pin,sA S 



(C4) 



In the all-pass configuration, the finesse of the cavity is 
given by 



F S = 



7T\/Pin,sAs 
1 - Pin,sA S ' 



(C5) 



Appendix D: Numerical cavity-reflectivity 
calculation 



For i G {5, P}, if Ti and Bi denote the forward and 
backward propagating fields, the intracavity field E{ can 
be written 



Ei(z)=Ti(z) + Bi(z). 



(Dl) 



With fi and bi, the slowly varying envelope amplitudes 
of the forward and backward propagating fields shown in 
Fig. [2^), we obtain 



Ei(z) = fi(z)e~ j ^ z + bi(z)e j ^ z , 



(D2) 



with ^ = 27rrid/\i. The field amplitudes are normalized 
in order to have I%{z) = \Ei(z)\ 2 . The calculation of the 
cavity reflection is a two point boundary value problem. 
It can be solved by a shooting method. The first bound- 
ary condition is that there is no incoming field from the 
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z > 0. This can be written by the following relation be- 
tween the forward and backward propagating field values 
at the back mirror 



bi(L) = pback,i/i(^)e 



-2j(3iL 



(D3) 



From this starting values we can deduce the values of the 
envelope amplitudes at the input mirror by integrating 
the following differential coupled equations 



dfp_ 

dz 
dbp 

dz 
dz 
dz 



~\ {(?p[ni(z) + 712(2)] + a p} fp( z ) 
- {(T P [m(z) + n 2 (z)} + a P } bp(z) 
-^KKW ~n 5 (z)} +a s }fs(z) 
7; Ws[ne(z) - n 5 (z)} + a s } b s (z), 



(D4) 



where the values of the NV center density are deduced 
from Eq. (|Alj) . We can obtain the input io,i = |^o,i| 2 
and output I out i — \E outji \ 2 intensities from 



— [/.«.)- «,A(o,] 



where K [n ^ for z G {P, S} («? n ^ + pf n i = 1) are the am- 
plitude mirror transmission coefficients. The calculation 
method consists in numerically optimizing the values of 
fi(L) to obtain the target values of Jo,i- The value of 
Rs = Iout,s/Io,s is then deduced with and without the 
microwave field applied. This is used to calculate the 
contrast C using Eq. (|8j) and the effective detected power 
max [i^s iis(0)] x P5. Finally, the minimum de- 
tectable magnetic field SB is evaluated using Eq. ([9]). 



TABLE I. Photophysical parameters of the six- level system 
sketched in Fig. [lji) . The transition rates kij are obtained by 
averaging data given in Ref. [32[. 1/T is the lifetime of level 
1 5). 71R is the spectral width of the 1042 nm zero-phonon line 
at room temperature. 



Parameter 


Value 


Reference 


Ap 


532 nm 


[19] 


Xs 


1042 nm 


[19] 


crp 


3 x 10" 21 m 2 


[33] 


/C31 — /C42 


(66 ± 5) fis- 1 


[32] 


&35 


(7.9 ±4.1) /is" 1 


[32] 


&45 


(53 ± 7) fis' 1 


[32] 


k 6 i 


(1.0 ±0.8) fis- 1 


[32] 


k 6 2 


(0.7 ±0.5) /is" 1 


[32] 


r 


1 ns" 1 


[18] 


7m 


« 2tt x 4 THz 





TABLE II. Physical parameters used in the single-pass NV - 
center IR absorption measurements 19 at room temperature. 
Optical losses are estimated from the transmission spectrum 
given in Ref. [34 ]. 



Parameter 


Value 


Reference 


n 


28 x 10 23 m~ 3 




T 2 * 


150 ns 


mi 


Ti 


2.9 ms 


[2S] (T = 300 K) 




400 MW/m 2 


[19] 


Iq,s 


10 MW/m 2 




Ps 


16 mW 




L 


300 /xm 


[19] 


n R 


2tt x 1.5 MHz 


[19] 


c 


0.003 


[19J (T = 300 K) 


as 


0.1-0.5 cm" 1 


[34] 



Appendix E: Tables 

The large value of T x shows that spin relaxation is 
negligible. Thus this is not taken into account in the 
rate-equation modeling of the six-level system. 



Appendix F: Sensitivity fundamental limit 

In this Appendix we derive the fundamental limit of 
the minimal detectable magnetic field value for methods 
based on IR absorption or visible fluorescence monitoring 
considering that the methods are limited by the photon 
shot-noise. 



1. IR absorption based magnetometer 

Using the expression of the ESR FWHM and Eq. (|2j) 
we obtain the following relation for the minimal de- 



tectable magnetic field 
SB = - 



(Fl) 



where iVph = Ps\s/(hc) is the number of detected IR 
photons per second. With Ns the number of IR photons 
collected per T| we have 



SB 



(F2) 



Now we estimate the maximal Ns value. Assuming an 
optimal contrast C = 1. When microwaves are switched- 
on, every photon is absorbed. We assume that one NV - 
center absorbs Ms IR photons per T 2 *. In many high- 
density samples, T 2 * < l/(&6i + ^62), and therefore we 
can consider that Ms < TT 2 *. We can thus write 



N s = M s ( 



(F3) 
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where N^ g is the number of NV centers in the singlet 
state when the microwave are switched-on and N^ g the 
number of NV - centers in the singlet for switched-off 
microwaves. This gives the number of photons which 
can be detected when the microwaves are switched-off 



N s = M S N 



i xP 45 



P35 



(F4) 



where N = nV is the number of centers with P35 = 
^35/(^35 + ksi) being the probability that NV - centers 
in level |3) (m s = 0) decay to the singlet and P45 = 
^45/(^45 + ^42) the probability that NV - centers in level 
1 4) (m s = ±1) decay to the singlet. The \ and | allow 
to take into account that only one quarter of the NV - 
centers are resonant with the microwaves^. We then 
have N s = K S M S N with 



^ _\( ^45 &35 

4 I &45 + &42 &35 + &31 



(F5) 



which is an approximated value for Rs(Qr) defined 
in section IIII Bl Note that if the IR power is such as 
IZs Ms > 1 the sensitivity is limited by the spin- noise. 



2. Fluorescence measurement based magnetometer 

For a magnetometer using the fluorescence signal mon- 
itoring and assuming that the ESR FWHM is 2/T|, the 
sensitivity is given b y 22 i 23 



SB { = 



(F6) 



where Cf is the contrast of the fluorescence signal and 
iVf the number of collected photons per T 2 *. When the 
microwaves are switched-off, the fluorescence signal is 



proportional to P31 = ksi/(ksi + fes) the probability 
that NV - centers in level |3) decay immediately to level 
|1). When the microwaves are switched-on the fluores- 
cence signal is proportional to P42/4 + 3P31/4 where 
P42 = £42/ (^42 + ^45) is the probability that NV - cen- 
ters in level |4) decay to level |2). Assuming that P31 « 1 
(^35 <C £31), the contrast Cf is given by 



Cf= 4 



*42 



C31 



35 



^42 



k 



45 



(F7) 



The number of collected photons per T 2 * is iVf = rjNM? 
where rj is the collection efficiency and Mf the number 
of emitted photons per T 2 * by one NV - center. Since 
1/&35 < T£ we have M f < k 31 /k 3 5- 



3. Comparison 



The two techniques can be compared by calculating 



IKsMs 



5Bf ^ 

IB ^ C f V rjM £ 
where we assume C ~ 1. For /C35 <C k 3 i and fe^ 
we have IZs ~ Cf and thus 



7p 



1 



I Ms 

M f ' rjK s ' 



(F8) 

j ^457 

(F9) 



Note that with values recalled in Tab. UU we obtain IZs ~ 
8.5% which corresponds to the optimal case asuming a 
total spin polarization. We deduce that Ms < 11 and 
Mf < 8. Assuming that M5 = Mf and considering a 
high value of the collection efficiceny (77 w 0.47 has been 
reported in Ref. [HI) we obtain SBf/SB w 5. 
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